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Abstract 

Using the two-loop approximation of chiral perturbation theory, we calculate the mo- 
mentum and density dependent single particle potential of nucleons in isospin-symmetric 
nuclear matter. The contributions from one- and two-pion exchange diagrams give rise to 
a potential depth for a nucleon at rest of C/(0,A;jo) = — 53.2MeV at saturation density. 
The momentum dependence of the real part of the single particle potential U{p,kfQ) is 
non-monotonic and can be translated into a mean effective nucleon mass of M* ~ 0.8M. 
The imaginary part of the single particle potential W{p,kf) is generated to that order 
entirely by iterated one-pion exchange. The resulting half width of a nucleon hole-state at 
the bottom of the Fermi sea comes out as VF(0,fc/o) = 29.7 MeV. The basic theorems of 
Hugenholtz- Van-Hove and Luttinger are satisfied in our perturbative two-loop calculation 
of the nuclear mean field. 



PACS: 12.38.Bx, 21.65.+f 
Keywords: Effective field theory at finite density. Real and imaginary part of the single particle 
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1 Introduction and summary 

The shell model and the optical model have always played a central role in the description of 
nuclear structure and nuclear reactions. The knowledge of single particle properties in dense 
nuclear matter is therefore of basic importance for the understanding of both nuclear structure 
and nuclear dynamics. Quantities like the nuclear mean field (or equivalently the single particle 
potential), the nucleon mean free path, the effective mass, etc., enter in the description of static 
as well as dynamic properties of nuclei (for a review see ref . [|I| ) . Considerable progress has been 
made during the past decades in calculating these single particle quantities from first principles, 
based on non-relativistic Brueckner theory |] and relativistic many-body theory ^ . Despite 
the widely different mechanisms which build up the nuclear mean field in these two theoretical 
approaches, the depth of the single particle potential comes out at approximately the same value 
when both calculations are restricted to first order in the renormalized in-medium NN-interaction 
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(G- matrix). Other properties, such as the momentum dependence of the nuclear mean field are 
also qualitatively similar in those two approaches. 

In a recent work ||^, we have used chiral perturbation theory for a systematic treatment of 
the nuclear matter many-body problem. In this calculation the contributions to the energy per 
particle, E{kf), originate exclusively from one- and two-pion exchange between nucleons and they 
are ordered in powers of the Fermi momentum kf (modulo functions of kf/m.,^ with denoting 
the pion mass). It has been demonstrated in ref. that the empirical saturation point and the 
nuclear matter compressibility K ^ 250 MeV can be well reproduced at order 0{kj) in the chiral 
expansion with just one single momentum cut-off scale of A ~ 0.65 GeV which parametrizes all 
necessary short range dynamics. Most surprisingly, the prediction for the asymmetry energy, 
Aq = 33.8 MeV, is in very good agreement with its empirical value. Furthermore, as a nontrivial 
fact pure neutron matter is predicted to be unbound and the corresponding equation of state 
agrees roughly with that of sophisticated many body calculations for low neutron densities 
Pn < 0.25 fm-3. 

Given the fact that the bulk properties of nuclear matter can be well described by chiral pion- 
nucleon dynamics treated up to three loop order it is natural to consider in a next step quantities 
which characterize in more detail the behavior of a single nucleon moving in the dense nuclear 
matter environment. It is the purpose of this work to calculate, using the same framework as 
in ref. [^, the momentum and density dependent (complex- valued) single particle potential of 
nucleons in isospin-symmetric nuclear matter. We will present here analytical expressions for 
both the real part U{p, kf) and the imaginary part W{p, kj) of the nuclear mean field generated 
by chiral one- and two-pion exchange. Our results can be summarized as follows: 

i) At equilibrium nuclear matter density chiral one- and two-pion exchange generate an at- 
tractive nuclear mean field for nucleons at rest, with a depth of t/(0, fc/o) = —53.2 MeV. 
This value is in very good agreement with the depth of the empirical optical model poten- 
tial Uq ~ —52 MeV deduced by extrapolation from elastic nucleon-nucleus scattering data 
d^, |[. In fact, the average nuclear potential employed in shell model calculations of nuclei 
has approximately the same depth [^]. 

ii) The momentum dependence of the real part of the single particle potential U{p,kfQ) is 
non-monotonic in the momentum interval < p < kfQ. Its maximal increase of AU = 
8.8 MeV can be translated into an average effective nucleon mass of M* = 0.82M. This is 
comparable to the empirical value of the effective nucleon mass in nuclear matter, M*^^ = 
(0.7 — 0.8)M, derived from experimental data in the framework of non-relativistic shell or 
optical models [|10], 0, g. 

iii) To the order considered here the imaginary part of the single particle potential W{p, kf) is 
entirely generated by iterated one-pion exchange. For a nucleon hole-state at the bottom of 
the Fermi sea, one finds a half width Thoie/'^ = W^(0, fc/o) = 29.7 MeV. This value is not far 
from results obtained in Brueckner-Hartree-Fock calculations based on realistic NN-forces, 
where W{Q, kfo) ~ 40 MeV 0. 

iv) The Hugenholtz- Van-Hove theorem [|I^] which relates the total single particle energy at 
the Fermi surface (i.e. the Fermi energy) to the nuclear matter equation of state is satisfied 
in our two-loop chiral perturbation theory calculation. Luttinger's theorem [jl4| does also 



hold since the calculated imaginary part, W{p,kf) ~ {kf — p)^, vanishes quadratically 
near the Fermi surface. The associated curvature coefficient C{kf) grows with density 
approximately as p^^^. 
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As a combined conclusion of ref . and this work one can thus state that perturbative chiral 
pion nucleon dynamics gives reahstic nuclear binding and saturation as well as realistic (in- 
medium) single particle properties. 



2 Single particle potential of nucleons in nuclear matter 

The nuclear mean field (or single particle potential) is generally defined through the momentum 
and density dependent nucleon self energy in (isospin symmetric) nuclear matter. In a diagram- 
matic calculation this quantity is derived from the in-medium nucleon propagator taking into 
account interactions among nucleons up to a certain order as well as Pauli-blo eking effects. On 
the other hand, the ground state energy density of nuclear matter is represented diagrammat- 
ically by analogous closed vacuum graphs which result from closing the nucleon line of a self 
energy diagram . As outlined in section 2 of ref.[P| the diagrammatic calculation of the ground 
state energy density can be organized in the number of so-called medium insertions. The latter 
is a technical notation for the difference between the vacuum and the bare in-medium nucleon 
propagator. Let us consider the structure of the energy density according to such an ordering 
scheme [Q. It consists of a sum of convolution integrals of the form, 

S[d] = J d^pi /Ci rf(pi) + J d^pid%IC2d{pi)d{p2) + J d^pid^p2d^p3}C3d{pi)d{p2)d{p3) 

+ J d^pid%d%d^p4 /C4 d{pi)d{p2)d{p3)d{p4) . (1) 

The one-body kernel JCi = 4Tk{\pi\) is four times the relativistically improved kinetic energy 
(see eq.(3)). /C2,3 are two- and three-body kernels related to contributions of closed diagrams 
with two and three medium insertions. In three-loop approximation the four-body kernel /C4 
is proportional to 6^{pi + P2 + P3 + Vi) and purely imaginary. The quantity d{pj) denotes 
the density of states in momentum space. Inserting the density of states of a filled Fermi sea, 
d{pj) = (27r)~^ 9{kf — \pj\), into eq.(l) one gets the energy density of nuclear matter, pE{kf), with 
the nucleon density p = 2A;j/37r^. The single particle potential can now be directly constructed 
from the energy density functional eq.(l) by adding a test nucleon of fixed momentum p to the 
filled Fermi sea. This situation is described by the density of states d{pj) = {2n)~^ 6{kf — \pj\) + 
T)S^{p — Pj) with the infinitesimal parameter 77 to be interpreted as (±) the inverse (infinite) 
volume. The plus sign applies for a particle {\p \ > kj) and the minus sign for a hole {\p \ < kf). 
Inserting this density of states into eq.(l) leads to 

S = pE{kj) + 4r] {Tk{p) + Uip, kf) + I W{p, kj)] , (2) 
with ^ ^ 

^'^^^ = h-m-^^ 

the relativistically improved kinetic energy. The factor 4 in eq.(2) simply counts the spin- and 
isospin multiplicity of a nucleon. The real and imaginary parts of the single particle potential 
have, according to eq.(l), a decomposition into two-, three- and four-body contributions, 

U{p, kf) = U2ip, kf) + U3ip, kf) , Wip, kf) = W2ip, kf) + W3ip, kf) + W,ip, kf) , (4) 

where the index on each term refers to the corresponding kernel /C„ in eq.(l). The energy per 
particle E{kf) in symmetric nuclear matter can also be reconstructed from these components of 
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the single particle potential by integrating over a Fermi- sphere of radius kj, 



kj 



dpp^ 



n{p) + lu2{p,kf) + ^u,{p,kf) 



(5) 



Note the weighting factors 1/2 and 1/3. These take care of combinatoric factors and the fact 
that three topologically equivalent closed vacuum diagrams with three medium insertions are 
generated by closing the nucleon line of a self energy diagram. Since this closing procedure 
must lead to a real energy per particle E{kf) the components of the imaginary part of the single 
particle potential have to fulfill the zero sum rule 



k} 



dpp^ 



\W2{P. kf) + ^W,ip, kf) + jW.ip, kf) 



0, 



(6) 



A further constraint on the single particle potential is given by the Hugenholtz- van- Hove theorem 
It states that the total single particle energy at the Fermi surface, p = kf, (i.e. the Fermi 
energy) is equal to the chemical potential. According to a general thermodynamic relation this 
chemical potential is the derivative of the energy density pE{kf) with respect to the particle 
density p. The Hugenholtz-van-Hove theorem reads explicitly, 



T,{kf)+U{kf,kf)=E{kf) + 



kf dEjkf) 
3 dkf 



W{kf, kf) = , 



(7) 



Note that the derivative term dE{kf) / dkf vanishes at nuclear matter saturation density and 
therefore the Fermi energy of equilibrated nuclear matter is equal to the negative average binding 
energy per particle, E{kfQ). Eqs. (5,6,7) serve as important numerical and analytical checks on 
our chiral perturbation theory calculation to be presented in the next two sections. 



3 Real part of the single particle potential 




Fig. 1: One-pion exchange Fock- diagram and iterated one-pion exchange Hartree- and Fock- 
diagrams. The combinatoric factor of these diagrams are 1/2, 1/4 and 1/4, in the order shown. 
Each nucleon propagator consists of a vacuum part and a medium insertion. 

In this section we present analytical results for the real part of the single particle potential 
U{p, kf) as given by chiral one- and two-pion exchange. The closed vacuum diagrams related to 
one-pion exchange (Fock- diagram) and iterated one-pion exchange (Hartree- and Fock-diagrams) 
are shown in Fig. 1. Each nucleon propagator consists of a vacuum part and a medium insertion 
(for details see eq.(3) in ref. 0). Self energy diagrams are obtained from these closed graphs by 
opening one nucleon line with a medium insertion. In order to keep their number small and in 
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order to avoid repetitions we associate the contributions to the (real) nuclear mean field U{p, kf) 
with the closed vacuum diagram before opening of a nucleon line. The index n on f/„(p, kf) (see 
eq.(4)) then indicates the number of medium insertions of the closed vacuum diagram before 
line opening. From a technical point of view, the calculation of a contribution Un{p,kf) (or 
Wn{p, kf)) involves an integral over the product of n — 1 Fermi spheres of radius kf. 

We start with the Ivr-exchange Fock-diagram in Fig. 1 with two medium insertions. With 
inclusion of the relativistic correction one finds a density and momentum dependent real 

two-body potential of the form, 
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where we have introduced the abbreviations u = kf/mT^ and x = p/m.,^. As in ref.Q, we choose 
the value qa = 1-3 for the nucleon axial vector coupling constant, f-^ = 92.4 MeV denotes the 
weak pion decay constant and = 135 MeV stands for the (neutral) pion mass. The iterated 
Ivr-exchange Hartree-diagram (second graph in Fig. 1) with two medium insertions gives rise to 
a real part of the single particle potential which reads 
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In this expression we have omitted the contribution of a linear divergence proportional to the 
cut-off A (see ref. [^). All such powerlike terms in A are collected in eq.(17). The iterated Ivr- 
exchange Fock-diagram (third graph in Fig. 1) with two medium insertions gives rise to a further 
contribution to the two-body potential U2{p, kf) of the form. 
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where we have again transferred a term linear in the cut-off A to eq.(17). The notation in eq.(lO) 
is to be understood such that the factor in the second line belongs to both integrals. Next, we 
consider the iterated Ivr-exchange Hartree-diagram (second graph in Fig. 1) with three medium 
insertions. In this case one gets three different contributions to the three-body potential f/3(p, kf) 
corresponding to the three possibilities of opening a nucleon line. Altogether they read. 
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with the auxihary functions, 



xy 



(12) 



Note that the expressions in eq.(ll) originate from a six-dimensional principal value integral 
over the product of two Fermi spheres of radius kf. The actual integrands have simple poles 
located on a five- dimensional quadratic hypersurface. The same features hold for the iterated 
iTT-exchange Fock-diagram (third graph in Fig. 1) with three medium insertions. Its contribution 
to the three-body potential U^i^p, kf) reads 
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where we have introduced some new auxiliary functions, 
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(14) 
(15) 



;i+x2-e2)2 + 4^2(l_^2)_ 

Note that all terms from iterated Ivr-exchange are proportional to the (large) nucleon mass, 
M = 939 MeV. The additional diagrams of irreducible 27r-exchange (not shown here, but see 
Fig. 4 in ref. pj) with two medium insertions lead to a real single particle potential which can be 



written as, 
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In order to avoid very lengthy analytical expressions we have stayed at a one-parameter integral 
representation in eq.(16). 

Finally, we give the complete expression for the power divergences specific to cut-off regular- 
ization, 



f/2(p, A;, 
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The term linear in the cut-off A stems from iterated Ivr-exchange with a contribution of the 
Hartree- and Fock diagram in the ratio 4:1. The term quadratic in the cut-off, on the other hand, 
originates from irreducible 27r-exchange. Note that the momentum independent contribution to 
the two-body potential f/2(p, kf) in eq.(17) is just twice its contribution to the energy per particle 
E{kf). This relative factor of 2 is typical for a momentum independent NN-contact interaction, 
to which the power divergences are completely equivalent, as emphasized in ref. H. 
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3.1 Results 



For the numerical evaluation of the real single particle potential U{p,kf) we use consistently 
the same parameters as in our previous work 0. There, the cut-off scale A = 646.3 Me V has 
been fine-tuned to the binding energy per particle, —E{kfQ) = 15.26 MeV. With this input the 
prediction for nuclear matter saturation density was po = 0.178 fm~'^ (corresponding to a Fermi 
momentum of kfo = 272.7 MeV) and the nuclear compressibility came out as = 255 MeV. 

In Fig. 2, we show the total real single particle potential U{0,kf) of our calculation for a 
nucleon at rest {p = 0) as a function of the nucleon density p = 2A;j/37r^. The shape of the curve 
in Fig. 2 is very similar to the nuclear matter equation of state (the so-called saturation curve). 
In comparison to the energy per particle, E{kf), the scale on the ordinate is stretched by a 
factor of about 3.5. Interestingly, the potential depth U{0, kf) reaches its minimum close to the 
saturation density po = 0.178 fm~^. The actual value at that point is f/(0, fc/o) = —53.2 MeV. 
This value is in good agreement with the depth of the empirical optical model potential Uq — 
—52 MeV deduced by extrapolation from elastic nucleon-nucleus scattering data 0, ||. The 
average nuclear potential used as input in shell model calculations of nuclei has approximately 
the same depth (see section 2.4 in ref.p|). For comparison, the calculation of ref.0 based on 
the phenomenological Paris NN-potential finds a potential depth of U{0,kfQ) ~ —64 MeV. In 
the relativistic Dirac-Brueckner approach of ref.p using the Bonn-A NN-potential a somewhat 
deeper real single particle potential with U{0,kfo) ~ —80 MeV has been found. It should also 
be noted that the potential depth obtained in the present work results from a cancelation of 
individually large attractive and repulsive terms. For example, when ordering the contributions 
to U{0, kfo) in third, fourth and fifth power of small momenta one has U{0, kfo) = (—314.7 -|- 
277.4 — 15.9) MeV. The decomposition of the same number into contributions from Ivr-exchange, 
iterated Ivr-exchange and irreducible 27r-exchange reads U{0, kfo) = (33.0 — 109.3 + 23.1) MeV. 
Similar features hold also for the binding energy per particle —E{kfo) = 15.26 MeV (see the 
discussion in section 2.5 of ref. 0). 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
P [fm '] 



Fig. 2: The real part of the single particle potential U{0, kf) at nucleon momentum p = versus 
the density p = 2A;^/37r^. 
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We note that both the sum rule eq.(5) and the Hugenholtz- van- Hove theorem eq.(7) hold 
with very high numerical accuracy in our calculation. The validity of the Hugenholtz-van-Hove 
theorem eq.(7) can easily be proven in the present framework on the basis of the density of states 
d{pj) = {2Tr)~^6{kf — \pj\) +'r]6^{p — pj) with \p \ = kf. 

In Fig. 3, the solid line shows the momentum dependence of the real single particle potential 
U{p, kfo) at saturation density for momenta from zero up to the Fermi surface, < p < kfo = 
272.7 MeV. The dashed line in Fig. 3 represents the total single particle energy, Tk{p) + U{p, kfo), 
i.e. the sum of single nucleon kinetic and potential energy. As required by the Hugenholtz-van- 
Hove theorem the dashed line ends at the Fermi surface p = fc/o = 272.7 MeV with the value 
E{kfo) = —15.26 MeV. The total single particle energy Tfc(p) + U{p,kfo) rises monotonically 
with the nucleon momentum p, as it should. Note however that there is no a priori guarantee 
for such a behavior in a perturbative calculation. The Hugenholtz-van-Hove theorem eq.(7) 
determines also the end point of the solid line in Fig. 3 as U{kfQ, kfo) = —54.0 MeV. Since 
it starts with almost the same value f/(0, fc/o) = —53.2 MeV at p = the real single particle 
potential U{p, kfo) cannot be a monotonically rising function of p in our calculation. As a matter 
of fact, one observes a downward bending of the real single particle potential U {p, kfo) for nucleon 
momenta above p = 180 MeV. Such a behavior is also found in the calculation of ref. based 
on the phenomenological Paris NN-potential, however, in much weaker form. The relativistic 
Dirac-Brueckner calculations of ref.|^ give, in contrast to this, single particle potentials which 
rise monotonically with the nucleon momentum p. 
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Fig. 3: The real part of the single particle potential U {p, kfo) versus the nucleon momentum p at 
nuclear saturation density, kfo = 272.7 MeV (solid line). The dashed line includes in addition 
the (relativistically improved) kinetic energy Tk{p) = p^/2M — p^/8M^. 



The momentum dependence of U{p,kfQ) translates into a variable effective nucleon mass 
M*{p) (the product of "/c-mass" and "i?-mass" divided by the free nucleon mass, according to 
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the nomenclature of ref.[^) by the relation, 

1 _ 1 1 dUjp, kfo) 
M*{p) M p dp ' ^ ' 

where M = 939 MeV is the (free) nucleon mass. Various many-body calculations |l|, find 
that the variable effective nucleon mass M*{p) is reduced below p ~ 0.8/c/o and enhanced around 
the Fermi surface p = kfo, typically by 20% on either side. The up- and downward bending of our 
real single particle potential U{p, kfo) (solid line in Fig. 3) also produces such an effect, but the 
momentum dependence is too strong. On the other hand, according to the sum rule eq.(5), the 
detailed momentum dependence of the components of the real single particle potential f/2,3(p, kj) 
and therefore also that of the variable effective nucleon mass M*{p) is of almost no relevance 
for the nuclear matter equation of state. A more suitable quantity for comparison is the mean 
effective nucleon mass, 

^ = ^ + ^ (19) 
M* k% ' ^ ^ 

with AU the maximal change of the potential U{p, kjo) in the momentum interval < p < kfQ. 
In the simple case that the real single particle potential grows quadratically with p, this M* agrees 
with the standard definition of the effective mass, eq.(18). By inspection of Fig. 3 one finds a 
maximal increase of our single particle potential U{p,kfo) of AU = 8.8 MeV. Via eq.(19) this 
shift AU translates into a mean effective nucleon mass of M* = 0.82M . This value is compatible 
with empirical effective nucleon mass M*^^ = (0.7 — 0.8)M derived from experimental data in 
the framework of non-relativistic shell or optical models |jl0|, |1T], 



The real single particle potential U{p,kf) has of course a smooth continuation above the 
Fermi surface, p > kf. We will not consider U{p,kf) in this region, for two reasons. First, 
most of the integral representations for the contributions U2,3{p,kf) given in section 3 do not 
hold anymore for p > kj. Secondly, when going above the Fermi surface the relevant nucleon 
momentum soon becomes too large for the application of chiral perturbation theory (which 
should be restricted to p < 300 MeV), keeping in mind that this theory is an expansion in small 
external momenta such as p and kf. 



4 Imaginary part of the single particle potential 

In this section, we discuss the imaginary part W{p,kf) of the single particle potential. For 
the reasons just mentioned, we will restrict ourselves to the momentum region below the Fermi 
surface, < p < kf. In this region nucleon-hole states with restricted momentum, < p < kf, 
can be prepared on top of the filled Fermi sea. Their total energy is —Tk{p) — U{p, kf)—iW{p,kf), 
according to eq.(2). The (positive) imaginary single particle potential W{p, kf) accounts for the 
finite life time of a nucleon-hole state via r^^^g = 2W{p, kf). By on-shell NN-scattering processes 
the energy of a deeply bound hole-state gets redistributed among two hole-states closer to the 
Fermi surface and a nucleon state in the continuum (i.e. above the Fermi surface). Energy and 
momentum conservation as well as Pauli blocking limit the available phase space which vanishes 
for p ^ kf as {kf — p)"^ according to Luttinger's theorem [Q. Within the present three-loop 
chiral perturbation theory calculation of nuclear matter the contributions to the imaginary single 
particle potential W{p,kf) arise entirely from iterated one-pion exchange. As done in section 
3 for the real part, we associate the contributions to W{p, kf) with the closed vacuum diagram 
(see Fig. 1) before opening of a nucleon line. Without going into further technical details, we 
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enumerate now the contributions to the components Wn{p,kf) and just indicate the type of 
iterated l7r-exchange diagram (Hartree or Fock) and the number of medium insertions n. 

i) Hartree-diagram with two medium insertions: 

+ (9 + + — - 2x'^) ln[l + {u + xf] + (9 + - — - 2x'^) 

X ln[l + {u- x)"^] + -(7 + Idu^ - ISx^) [ arctan(ii + x) - arctan(u - x) |(,20) 

X J 

with the abbreviations u = k^/rriT^ and x = plm-^. 

ii) Fock-diagram with two medium insertions: 



W^2(P,A;/) 



(47rA)4 I 4 ^ 

{u-x)l2 



X 

6 60 



Jo J(u-x)/2 IX 



1 + 2^2 



ln(l + 4a . (21) 



iii) Hartree diagram with three medium insertions: 



(47r/.)4 
h lOxy 



,2„.2n 



dy\ {u -X y 



arctan s — s + 



2s^ + 

1 + S2 



-21n(l + s2 



s I [y — 1) 



+ 



« 2^2 



X 



2(7^ + a' 



l + (72 



+ (3 + 2ii^ - 2x'y') - ln(l + s^) 



21n(l + (7^) 



(22) 



The auxihary functions s and cr have been defined in eq.(12). 
iv) Fock diagram with three medium insertions: 



dy 



1 4V1 -y'^-z' 



- ln(l + s^) ln(l + e) - r 



(23) 



The auxihary functions t and R have been defined in eq.(15). 
v) Hartree diagram with four medium insertions: 



QTig\Mm\ f 7 8x 
(47r/.)^ 



--- + --31n(l + 4.^) + (5a:^- 



15 



Ax 



arctan 2x 



+ 



dy 



IQxy 



— — arctan s 



+ (3 + 2u'^ - 2x'^y^)[ln{l + s^) - 



3^4 ^4^2(1 _ y2^ 



^ 2 ^ l + s2 



(24) 



vi) Fock diagram with four medium insertions: 



X 



- ln(l + a^) + 9{x -0(^1' + 



(25) 
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with a new auxiliary function 



w 



(26) 



The total sum of these six contributions eqs. (20-25) evaluated at zero nucleon momentum 
{p = 0) can even be written as a closed form expression, 



W{0,kf) 



37ig\Mmt 
(47r/.)4 



u 



+ 



6 + 



2 
ln(2 



u 



^ H h Li2(-1 - 

l + u2 12 ^ ' 



(27) 



with u = kf/m^ and Li2(— a^^) = Jq d( {C + 'j^)"^1ii(^ denotes the conventional dilogarithmic 
function. Note that eq.(27) determines the half width of a nucleon-hole state at the bottom of 
the Fermi sea. The density dependence of W{0, kj) in eq.(27) goes approximately as p^^^ since 
the first term, m^/2, in the curly bracket dominates. This is in fact also the only term which 
survives in the chiral limit rriTr = 0. 



4.1 Results 

Evidently, the results for the imaginary single particle potential W{0, kf) are completely parame- 
terfree (to the order we are working here). At saturation density po = 0.178 fm~^ (corresponding 
to kfQ = 272. 7 Me V) we find from eq.(27) a half width for a nucleon-hole state at the bottom 
of the Fermi sea of W{0,kf) = 29.7 MeV. This is not far from the value W{0,kf) ~ 40MeV 
obtained in the self-consistent Brueckner calculation of ref.[^ using the phenomenological Paris 
NN-potential. On the other hand the calculation of ref. employing the Gogny Dl effective 
interaction finds about half of the value of ref.|0], namely W{0,kf) ~ 20 MeV. Taking the re- 
sults of these two calculations as a reasonable benchmark our prediction for the imaginary single 
particle potential at zero-momentum, W{0, kf) = 29.7 MeV, can be considered as fairly realistic. 




100 150 
p [MeV] 



Fig. 4: The imaginary part of the single particle potential W{p, kfo) versus the nucleon momen- 
tum p at nuclear matter saturation density, kfQ = 272.7 MeV. 
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In Fig. 4, we show the momentum dependence of the imaginary part of the single particle 
potential W{p, kfo) at nuclear matter saturation density, kjQ = 272.7 MeV. As a consequence of 
the decreasing phase space available for the redistribution of a nucleon-hole state's energy the 
curve in Fig. 4 drops with momentum p and W{p, kfo) reaches zero at the Fermi surface p = kjQ. 
From the formulas eq. (20-26) given above one can actually proof that both the value W{kf, kf) 
and the slope [dW{p,kf)/dp]p=kj: vanish identically at the Fermi surface p = kf. This is even 
separately true for the contributions of the iterated Ivr-exchange Hartree- and Fock-diagram. 
Luttinger's theorem is therefore satisfied in our perturbative calculation. The imaginary 
single particle potential grows quadratically near the Fermi surface, 



wip,kj) = cikf)-ip~kjY + 



{2i 



with C{kf) a density dependent curvature coefficient. The quantity C{kf) is useful in order to 
characterize the damping of the single particle motion of valence nucleons in nuclei [0]. In Fig. 5, 
we show the dependence of the curvature coefficient C{kf) (given in units of GeV^^) on the 
density for p < 0.4fm~'^. One observes an approximate p^^^-dependence of C{kf). This behavior 
can be easily understood when taking the chiral limit m^r = 0. In that case the proportionality 
C{kf) ~ k'j becomes a simple consequence of (mass) dimension counting. 

Finally, we note that the zero sum rule eq.(6) for the two-, three- and four-body components of 
the imaginary single particle potential W^2,3,4(p, kf) holds with high numerical accuracy separately 
for the iterated Ivr-exchange Hartree- and Fock-diagram. This serves as an important check on 
our calculation. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
P [fm"] 



Fig. 5: The curvature coefficient C{kf) of the imaginary part of the single particle potential 
W{p^kf) versus the nucleon density p = 2kj^/37i^. 

5 Concluding Remarks 

In this paper, we have continued the systematic treatment of the nuclear matter many-body 
problem in the framework of chiral perturbation theory by calculating the (complex-valued) 
single particle potential or nuclear mean field. From the good results of our previous paper 
and the present work one can indeed conclude that perturbative chiral pion-nucleon dynamics 
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is able to produce nuclear binding and saturation (including its isospin dependence) as well as 
realistic in-medium single particle properties (including the average nuclear mean field). Most 
surprisingly, all this can be achieved with only one adjustable parameter, namely a single cut-off 
scale A ~ 0.65 GeV which represents short distance NN-dynamics. 

Our approach to the nuclear matter problem is in many respects different from most other 
commonly used ones. First, we do not start from a so-called realistic NN-potential (which fits 
deuteron properties and NN-phase shifts), but we use instead well-founded chiral pion-nucleon 
dynamics. From that we calculate analytically in-medium multi-loop diagrams up to a certain 
order in a small momentum expansion. This "inward-bound" method starts from large distances 
(small Fermi momentum kf) and systematically generates the pion-induced correlations between 
nucleons as they develop with decreasing distance (increasing kf) down to the average NN- 
distance in nuclear matter, d ~ 1.8 fm, which is just slightly larger than the pion Compton wave 
length = 1.46 fm. 

Of course, open questions concerning the role of yet higher orders in the chiral expansion 
remain and should be explored. 
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